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INTRODUCTION 


In  recent  years,  an  explosion  in  the  capabilities  o-f 
Computational  Fluid  Dynamics  (CFD)  has  occured.  Since  the 
introduction  o-f  the  Small  Perturbation  method  of  Murman  and  Cole, 
transonic  CFD  has  quickly  evolved  towards  Full  Potential  and  Euler 
methods  capable  of  analyzing  a  variety  of  three-dimensional 
configurations.  Although  much  Success  has  been  realized  in  this 
field,  continued  success  is  inevitable  with  the  ever  increasing 
advances  in  computing  power  and  faster  algorithms.  The  purpose  of 
this  research  is  to  study  the  applicability  of  pseudo-spectral 
schemes  for  transonic  CFD. 

MODEL  PROBLEMS 

Several  time-stepping  procedures  were  investigated  to  become 
familiar  with  techniques  appropriate  for  the  pseudo-spectral  scheme. 
The  procedures  were  tested  on  two  model  problems  (including  Burgers' 
inviscid  equation),  and  each  procedure's  stability  criterion  was 
determined.  The  main  purpose  of  these  exercises  was  to  provide  the 
author  experience  and  a  sense  of  confidence  in  several  time-stepping 
procedures,  specifically,  their  accuracy,  and  their  stability  as 
defined  by  the  mathematical  analyses. 

The  first  model  problem  was  to  solve  a  simple  differential 
equation  with  a  known  solution  for  comparison.  The  differential 
equation  simulated  was 


=  1  ,  ( 1 ) 


dy/dt  -  a*y 


exp(b*t>,  with  initial  condition  y(0) 


where  y  is  a  -function  o-f  t,  and  a  S<  b  are  complex,  constants.  The 
exact  solution  to  equation  (1)  is 
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is 


y ( t )  =  ( 1/ (b-a) ) *exp <b*t>  +  (1  -  1 / (b-a ) ) *exp < a*t ) . 


The  tour  time-stepping  procedures  used  were  Euler's  method,  a 


backward  method,  the  mid-point  rule,  and  a  leap-frog  method.  These 


methods  and  their  stability  are  discussed  below. 


Euler's  method  is  a  f orward-dif f erence  scheme  which  is 


first-order  accurate  in  time.  This  scheme  is  described  as 


Y<n+1)  =  (  1  +  a*d)*Y(n)  +  d*exp  (b*n*d )  , 


where  Y  is  the  approximation  of  y,  d  is  the  time  step,  and  n  is  the 


time  level,  ie.  t  =  n*d.  A  stability  analysis  of  this  scheme  defines 


the  boundary  on  the  time  step  as 


d  <=  -2*a 1 / ( a 1 **2  +  a2**2) , 


where  a  =  al  +  i*a2,  and  i  is  the  square  root  of  -1.  Thus,  from 


luation  (4)  it  is  easily  seen  that  Euler's  method  cannot  be  used  if 


al  is  greater  than  or  equal  to  0,  or  if  a2**2  >>  absCall 


The  backward  method  is  a  backward-di f f erence  scheme  which  is 


also  first— order  accurate  in  time  and  is  described  as 


(1  -  a*d ) *Y (n  +  1 )  =  Y(n)  +  d*exp < b* (n+ 1 ) *d )  . 


A  stability  analysis  indicates  that  this  scheme  becomes  more  and  more 
stable  as  the  time  step  goes  to  infinity.  However,  since  this  scheme 
is  only  -first-order  accurate  in  time,  only  small  time  steps  should  be 
appl i ed . 

The  mid-point  rule  is  a  second-order  time  accurate  scheme 
described  by 

(1  -  .  5*a*d>  *Y  (n  +  i  )  =  (1  +  .5*a*d)*Y<n)  +  d*exp  (b*  <n+.  5>  *d )  .  (6) 

The  stability  of  this  scheme  requires  al  to  be  strictly  less  than  0. 

The  leap-frog  method  is  also  second-order  accurate  in  time  and 
is  described  as 

Y (n  +  1 )  *  Y(n-l)  +  2*a*d*Y(n>  +  2*d*exp <b*n*d) .  (7) 

Starting  the  leap-frog  scheme  also  requires  an  additional  unknown, 
ie.  Y (1 ) ,  which  can  be  approximated  using  the  mid-point  rule  so  that 
second-order  accuracy  is  maintained.  The  leap-frog  scheme  is  stable 
i  -f  d  <  1/al,  thus  useless  i -f  al  is  negative. 

The  above  time-stepping  procedures  were  tested  -for  a  variety 
of  complex  constants  a  &  b  to  verify  the  stability  estimates  and 
order  of  accuracy  determined  by  the  mathematical  analyses.  With  the 
confidence  provided  by  the  first  simple  model  problem,  it  was  time  to 
study  an  equation  with  the  flavor  of  a  transonic  flow  problem. 

The  second  model  problem  investigated  was  Burgers'  inviscid 
equation.  In  conservation  form,  this  quasi-linear  partial 


di  i  -f  erenti  al  equation  is  described  by 


du/dt  +  . 5*d (u**2) /dx  =  0, 

i 


or  in  quasi-linear  -form,  described  by 


du/dt  +  u*du/dx  =  0. 


(8a) 


(8b) 


Let  us  consider  the  initial  conditions 


U  (X 

, 0 )  =  0.5, 

x  <  0, 

U  (X 

,0)  =  1.25  -  sqrt (0. 5625  +  2*x  -  x**2) , 

0  < 

:  x  <  i. 

(8c) 

u  (x 

,  0 )  =  0 .  0 , 

1  < 

x  . 

The  initial  conditions  were  chosen  such  that  a  "shock"  would 
eventually  -form  in  the  solution,  thus  providing  a  reasonable  test 
case  -for  a  scheme  to  be  used  -for  transonic  -flow  calculations.  The 
exact  solution  to  equations  (8)  is  presented  and  discussed.  In 
addition,  several  numerical  methods  are  used  to  calculate 
approximations  -for  the  above  problem  and  are  compared  to  the  exact 
solution  through  discussion.  The  numerical  methods  are  analyzed  -for 
stability  criteria  and  order  o-f  dissipation  when  applied  to  a  -form  o-f 
equations  (Bb)  with  constant  coe-f  -f  i  c  i  ents ,  ie.  du/dt  +  A*du/dx  =  o , 
where  u  is  a  vector  and  A  is  a  constant  matrix.  The  stability 
criteria  used  (-for  approx  i  mat  i  ng  the  di  t -f  erent  i  al  equation)  is  the 
most  restrictive  one  obtained  by  allowing  the  constant  A  to  vary  over 
the  interval  containing  all  values  o-f  u,  le.  CO.  0,0. 53. 


Observe  that  equation  (8b)  can  be  solved  using  the  method  of 
characteri  sties.  Thus,  u  =  u(r,s),  where  r  describes  the  x-location 
on  the  data  line  where  the  characteri stic  line  starts  and  s  is  in  the 


direction  o-f  the  characteri  sti  c  line.  Note  that  at  some  time,  the 
characteristic  lines  begin  to  merge  to  -form  a  shock.  The  movement  of 
this  discontinuity  can  be  described  by  the  entropy  jump  condition 
(see  Reference  1).  By  observation  of  equation  (8b)  and  initial 
conditions  (8c) ,  it  can  be  noted  that  a  shock  will  surely  farm  and 
will  travel  to  the  right.  The  presentation  of  the  exact  solution 
will  be  given  in  three  parts,  specifically,  the  time  before  the  shock 
(ie.  0  <  t  <  tl),  during  the  growth  of  the  shock  (ie.  tl  <  t  <  t2) , 
and  after  the  shock  is  fully  mature  (ie.  t2  <  t).  Using  the  method 
of  characteristics,  the  exact  solution  for  the  first  time  interval 
can  be  shown  to  be 


k 

> 

» 


u  (x 

,t) 

=  0.5, 

x  <  t/2, 

U  (  X 

,t) 

=  u  ( r  ,  0 )  , 

t/2  < 

x  <  1,  (9a) 

U  (  X 

,t) 

1! 

O 

■ 

o 

«• 

1  < 
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where  u(r,0)  is  defined  by  (8c)  and  r  is  given  by 

1 . 25*t  -  sqrt (0.5625  +  2*r  -  r**2)*t  +  r  -  x  =  0.  ■ ?b ) 

! 

Solutions  (9)  are  valid  only  on  the  time  interval:  0  <=  t  =  tl, 

where  tl  is  the  time  at  which  the  shock  first  forms.  The  value  of  tl 
can  be  determined  by  finding  the  first  occurrence  of  abs(du/dx> 
approaching  infinity,  ie.  when  character  i  stic  lines  begin  to  merge 


together.  These  calculations  show  that  the  shock  begins  at  x  =  0.375 
and  t  =  0.75,  thus  de-fining  tl. 

The  exact  solution  -for  the  second  interval  is  described  by 


U  (  X  ,  t  ) 

=  0.5, 

X  < 

x shock , 

U  (X  ,  t  ) 

=  u  ( r  ,  0 )  , 

x shock  <  x  < 

1, 

( 10) 

U  (  X  ,  t  ) 

=  0 .  o ,  • 

1  <  X, 

where  u(r,0>  and  r  are  as  before,  and  xshock  is  the  x-location  o-f  the 
shock.  To  de-fine  the  motion  o-f  the  shock  during  its  growth,  a 
numerical  time-stepping  procedure  was  utilised  based  on  conserving 
properties  across  the  discontinuity.  This  shock  motion  is  provided 
i n  Table  I  where  shock  location,  r — ,  r+,  u+,  and  shock  speed  are 
provided  as  a  -function  o-f  time.  (Here,  the  +  and  -  indicate  on  which 
side  of  the  shock  the  property  is  given.)  Note  that  the  shock  speed 
monoton  i  cal  1  y  decreases  -from  0.5  to  0.25  during  its  growth.  Another 
important  item  to  note  -from  Table  I  is  the  time  of  shock  maturity, 
which  defines  t2  to  be  2.7953-3. 

The  exact  solution  for  the  third  interval  is  now  known  to  be 
described  by 


u  ( x , t )  =  0.5, 


x shock  , 


u  ( x , t )  =  0.0, 


x shock 


(11) 


The  exact  solution  to  Burgers'  problem  (8)  is  described  by 


equations  (9-11)  and  is  illustrated  in  Table  I  and  Figures  1-2. 

Figure  1  provides  a  solution  profile  for  several  "key"  times.  The 
first  time  (t  =  0.0)  illustrates  the  initial  conditions,  the  second 
time  (t  =  tl  =  0.75)  indicates  the  start  of  the  shock,  the  third 
fourth  times  (t  =  1.5  ?<  2.25)  help  show  the  growth  of  the  shock,  and 
the  fifth  time  (t  =  2.79533)  shows  the  instant  the  shock,  becomes 
fully  mature.  Figure  2  illustrates  the  position  of  the  shock  as  a 
function  of  time.  Notice  that  the  shock  forms  at  x  =  0.375  and  t  = 
tl  =  0.75,  and  becomes  fully  mature  at  x  =  1.0  and  t  =  t2  =  2.79533. 
Also,  notice  that  the  solution  to  the  left  of  the  shock  and  the  r  =  0 
char acter i st i c  is  u  =  0.5,  the  solution  to  the  right  of  the  mature 
shock  and  the  r  =  1  characteristic  is  u  =  0.0,  and  the  solution  in 
between  is  determined  by  equations  (9b),  (8c),  and  Table  I.  With  the 

exact  solution  of  Burgers  '  equation  (8)  known,  several  numerical 
approaches  will  be  evaluated. 

The  first  numerical  approx i mat i on  to  equations  (8)  uses 
MacCormack's  method  (Reference  2).  This  scheme  is  second-order 
accurate,  conditionally  stable,  fourth-order  dissipative,  and  can  be 
split  into  two  steps;  a  predictor  step  and  a  corrector  step.  For  the 
conservation  form  (8a),  the  scheme  is  described  by 

Vtmp(i)  =  V ( i , n )  -  q*0 . 5* ( V < i + 1 , n ) **2  -  V(i,n)**2) 

V(i,n+1)  =  0.5*(V(i,n)  -  q*0 . 5* ( Vtmp ( i ) **2  -  Vtmp ( i - 1 ) **2 ) )  (12a) 


+  0 . 5*Vtmp ( l ) 


and  -for  the  quasi-linear  -form  <8b>,  the  scheme  is  described  by 

Vtmp(i)  =  V(i,n)  -  V(i,n)*q*(V(i+l,n)  -  V(i ,n) ) 

V ( i ,  n  + 1 )  =  0.5*(V(i,n>  -  Vtmp ( i > *q* ( Vtmp ( i )  -  Vtmp(i-l))  (12b) 

+  0. 5*Vtmp ( i ) , 

where  V  is  the  appro:;  i  mat  i  on  o-f  the  exact  solution  u,  Vtmp  is  a 
temporary  value  determined  by  the  predictor  step,  i  indicates  the 
spatial  location  (ie.  x  =  i*delta-x),  n  indicates  the  time  level  (ie, 
t  =  n*delta-t)  ,  and  q  is  the  ratio  o-f  delta-t  to  delta-x.  As 
mentioned  previously,  MacCormack's  scheme  is  conditionally  stable  and 
through  mathematical  analyses,  stability  requires  q  <  =  2.0.  Figure  3 
illustrates  the  solution  o-f  equations  (12a)  -for  initial  data  (8c)  at 
times  o-f  0,  1,  2,  3,  4,  5,  and  6.  Notice  that  the  shock  is  -fully 
mature  by  time  equal  to  3.0.  (Recall  that  the  exact  solution  has  the 
first  occurrence  of  a  mature  shock  at  time  equal  to  2.79533.)  Also, 
notice  that  the  mature  shock  speed  is  almost  exactly  equal  to  0.25, 
which  compares  very  nicely  with  the  exact  solution.  Another  item  to 
notice  in  Figure  3  is  that  the  solution  seems  to  be  stable  with  no 
pre-shock  oscillations  evident.  This  solution  was  obtained  with 
delta-x  =  0.05,  and  delta-t  =  0.10,  which  by  the  mathematical 
analysis  is  neutrally  stable  (ie.  q  =  2.0).  Figure  4  shows  the 
solution  of  equations  (12b)  for  the  same  time  levels.  This  solution 
is  stable;  however,  the  most  important  thing  to  note  is  that  the 
shock  stagnates  at  x  =  1.0.  This  figure  definitely  indicates  the 
importance  of  using  the  conservation  form  of  Burgers'  equation  as 
opposed  to  using  its  quasi— linear  form.  Another  item  to  notice  is 


that  a  pre-shock  oscillation  seems  to  be  -forming  as  the  time 


increases.  This  solution  also  used  delta-:-;  =  0.05  and  delta-t  = 

0.  10. 

The  second  approach  used  to  solve  Burgers'  problem  (8)  was  an 
implicit  scheme  developed  by  Beam  and  Warming  (Reference  3).  This 
scheme  has  been  analyzed  to  be  neutrally  stable  and  non-di ssi pati ve. 
For  convenience,  delta-;<  and  delta-t  have  been  chosen  to  be  equal  to 
that  used  -for  MacCormack ' s  method.  This  scheme  approximates  the 
conservative  -form  (3a)  by, 

0. 25*q* (V (i+1 ,n)*V(i  +  l ,n  +  l>  -  V ( i -1 , n ) *V ( i -1 ,n+l ) ) 

+  V(i ,n+l>  =  V ( i , n ) ,  (13a) 


and  the  quasi-1  inear  -farm  ( Qb )  by, 

V(i,n+1)  +  0.5*q*V(i ,n)*(V(i+l ,n  +  l)  -  V(i-l,n+l>)  =  V(i,n).  (13b) 

Figure  5  shows  the  solution  o-f  (13a)  for  times  equal  to  0.0,  0.5, 

1.0,  1.5,  and  2.0.  Notice  the  extreme  pre-shock  oscillations  for 

this  solution.  This  occurs  because  the  Beam-Warming  scheme  is 
non-di ssi pati ve.  The  solution  to  equation  (13b),  shown  in  Figure  6, 
also  has  the  extreme  pre-shock  oscillation.  To  remedy  this  problem, 
a  fourth— order  dissipative  term  is  added  to  the  right-hand  side  of 
equations  (13).  This  term  is  explicitly  defined  as, 

-0. 125*w* ( V ( i +2 , n )  +  V(i-2,n) 

-  4*  ( V  (  l  +  1  ,  n  )  +•  V  (  i  -  1  ,  n  )  >  +  6*0  ( l  ,  n  i  )  'll' 


v  > n  .  *  \> * . .  v  - v 
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where  w  determines  the  amount  of  dissipation  and  is  allowed  to  have 
the  values  0  <=  w  <=  1  according  to  a  linearized  stability  analysis. 
Figure  7  shows  the  solution  of  equation  (13a)  with  the  dissipation 
term  (14)  tor  times  o-f  0.0,  0.5,  1.0,  1.5,  2.0,  2.5,  and  3.0  with  w 
equal  to  0.75.  The  dissipation  term  had  very  little  effect  on  the 
pre-shock  oscillations  when  applied  to  the  conservative  form  (13a), 
yet  damped  the  oscillations  of  the  quasi-1 inear  form  (13b)  very 
nicely  as  is  shown  in  Figure  8.  Also,  notice  in  Figures  7-8  that  the 
dissipation  term  produced  a  small  oscillation  at  the  base  of  the 
shock . 

The  third  approach  used  to  solve  Burgers'  problem  (8)  was  a 
pseudo-spectral  method.  Since  the  solution  to  Burgers'  inviscid 
problem  with  the  given  initial  conditions  is  obviously  non-periodic, 
a  Chebyshev  polynomial  expansion  is  employed  to  evaluate  the  spatial 
derivatives  in  a  spectral  manner.  In  particular,  the  author 
investigated  the  space  interval  (-1,3)  for  times  between  0  and  3. 

The  spatial  interval  is  discreatized  with  I  cells,  thus  the  solution 
is  described  by  1  +  1  values  at  the  collocation  points.  These  points 
are  located  at, 

X(i)  =  S*cos <p i *i / I )  +  T  f or  i  =  0 , 1 , . . . , I ,  (15> 

where  S  is  a  scale  parameter,  and  T  is  a  translation  parameter.  For 
the  spatial  interval  investigated,  S  =  2,  and  T  =  1.  The  solution  at 
the  I + 1  points  is  described  by, 

U(i)  =  Sum  (  j  =0  ,  I  )  A  (  j  )  *cos  (pi  *i  *  j  /  I  )  for  i  =  0,1,...  ,1,  1  1  ■=  1 


where  the  A(j)  coef  f  i  ci  ents  can  be  evaluated  by  the  inverse 
transform, 

A ( j )  =  2/ < I*C  <  j ) )  *  Sum (i =0,1)  U ( i > *cos (pi * j*i  / 1 ) /C ( i  )  ,  (17 

where  C(j)  =  2  -far  j  =  0  or  I,  and  C(j)  =  1  otherwise.  Equations 
(16-17)  can  be  evaluated  by  Fast  Fourier  Transforms  (FFT)  i f  I  i s  a 
composite  integer.  By  differentiating  equation  (16),  one  obtains, 

UX <  i  )  ~  Sum (j =0,1)  B< _j  >  *cos  (p  i  *i  *  j  / 1 ) /S ,  (IS 

where  B(j)  can  be  computed  from  A(j>  by  using  the  recurrent 
rel ati ons , 

B  ( I )  =  0  , 

B(I-l)  =  2*I*A ( I ) ,  (19 

B ( j  )  =  <B ( j+2)  +  2*( j  +  l)*A( j  +  1) > /C( j) . 

The  time-stepping  algorithm  used  is  a  simple  predi ctor-corrector 
scheme,  and  is  described  by, 

Utmp(i)  =  U(i,n)  +  k*U ( i , n ) *UX ( i , n ) , 

U(i,n+1)  =  0.  5*  <  U  (  i  ,  n  )  k*Utmp  (  i  )  *UXtmp  (  i  )  +  Utmp(i)),  (20 

where  k  is  the  delta-time  specified.  Since  the  solution  to  equation 
(Si  is  discontinuous,  spurious  oscillations  are  introduced  in  the 
pseudo— spectral  scheme  such  that  the  scheme  inevitably  diverges.  To 


aid  convergence,  a  spectral  -filter  is  applied  to  the  A(j>  terms  to 
reduce  the  amplitude  o-f  the  high  -frequency  contributions.  This 
■filtering  scheme  is  given  as, 

A ( j )  =  A ( j )  for  j  =  0,1,. . . ,J0, 

A(j)  =  A ( j ) *exp (RHO* ( ( j— JO) / ( I  — JO) ) **4 >  otherwise,  (21) 

where  RHO  is  the  decay  parameter  and  JO  is  the  cut-o-f-f  -filter 
parameter . 

The  above  pseudo-spectral  method  was  tested  on  Burgers ' 
equation  (8)  in  both  conservative  and  quasi-1 inear  form.  Figur6  9 
illustrates  the  pseudo-spectral  appro:-:  i  mat  ion  of  Burgers'  equation 
(8a)  in  conservation  form  for  several  times.  Notice  the  growth  of 
the  pre-shock  oscillations  with  time.  Figure  10  illustrates  the 
solution  of  the  quasi-1 inear  equation  (8b)  for  the  same  times. 

Notice  that  the  pre-shock  oscillations  are  more  pronounced  and  the 
shock  movement  is  slower  than  the  conservative  solution.  In  fact, 
the  quasi-1 inear  shock  will  stagnate  at  X  =  1,  while  the  conservative 
shock  will  continue  to  move  to  the  right  with  a  speed  appro;:  i  matel  y 
equal  to  the  exact  solution.  This  behavior  between  the  conservative 
and  quasi-linear  solutions  has  been  observed  before  in  the  MacCormack 
and  Beam-Warming  schemes. 

In  addition  to  the  above  effort,  some  time  was  spent  studying 
Fast  Fourier  Transform  algorithms  (Reference  4)  (which  are  essential 
to  the  economical  competitiveness  of  the  spectral  schemes),  spectral 
filtering  schemes,  and  artificial  dissipation  terms  (References  5-8) 
''which  are  neccessary  for  shoe)  -capturing  when  using  the  spectral 


methods) . 


After  experimenting  with  these  shoe k -c aptur i ng  techniques 
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and  examining  similar  results  -from  other  researchers  (including  those 
in  -favor  of  the  pseudo-spectr al  scheme),  it  seems  that  the 
pseudo-spectral  method  is  very  well  suited  for  smooth  solutions. 
However,  if  the  solution  contains  a  discontinuity,  the  high-order 
accuracy  of  the  pseudo-spectr al  scheme  can  quickly  degenerate  to  as 
poor  as  first-order  accuracy;  easily  eliminating  any  advantage  the 
pseudo-spectral  scheme  may  have  over  the  simpler  second-order 
f i ni te-di f f erence  or  finite-element  schemes.  The  substantial  loss  in 
accuracy  of  a  spectral  shock-captur i ng  scheme  is  verified  by  Hussaini 
(Reference  9)  when  studying  the  time  development  of  shock  waves  in  a 
spiral  galaxy.  In  order  to  retain  the  high  accuracy  of  the 
pseudo-spectral  method  (for  solutions  with  discontinuities),  a 
shock-fitting  technique  could  be  employed.  The  solution 
discontinuity,  by  being  fitted,  could  be  represented  as  a  boundary, 
and  by  using  Chebyshev  polynomials,  one  could  handle  the  non-periodic 
boundary  conditions  (Reference  10  studies  a  similar  approach). 
Although  this  concept  would  retain  the  high-order  accuracy  desired, 
it  would  also  introduce  a  new  set  of  problems.  These  problems  incldde 
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1)  the  prohibitively  small  time  steps  allowed  by  the  stability  of 
schemes  using  the  Chebyshev  spacing,  2)  problems  handling  multiple 
shocks,  and  3)  problems  expanding  the  shock-fitting  scheme  to 
three-dimensional  solutions.  In  view  of  the  obstacles  outlined 
above,  the  author  feels  that  the  pseudo-spectral  scheme  is  not  well 
suited  for  "practical"  transonic  problems. 


The  proposed  research  is  to  develop  a  transonic  Euler  method 
based  on  a  -finite-element  scheme  similar  to  that  of  Jameson  et.  al  . 


The  intention  is  to  develop  a  flow  solver  un i versal  to  any  type 
configuration  regardless  of  dimension  or  topology,  and  initially 
develop  a  grid  generation  and  post-pr ocessor  package  for 
two-dimensional  airfoil  geometries.  Although  the  flow  solver  will  be 
very  general,  the  grid  generation  and  flow-solution  post-processor 
packages  will  be  particular  to  each  case,  thus  future  work  will 
include  expanding  the  capability  of  the  method  by  developing  grids 
for  multielement  airfoils  and  three-dimensional  configurations. 

Other  items  which  may  be  explored  are:  1)  multigrid,  2)  "smart" 
dissipation  at  shock  locations,  3)  code  designed  for  parallel 
processing,  4)  capability  for  an  inverse  mode,  5)  viscous  effects,  6) 
local  mesh  refinement,  7)  inclusion  of  power  effects,  and  8)  more 
general  grid  generation  tools  to  ease  future  geometric  needs. 


CONCLUSIONS 

*  The  pseudo— spectral  scheme  is  excellent  -for  smooth 
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* 
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solutions. 

The  pseudo-spectral  scheme  can  easily  degenerate  to 
f  i  rst-or  der  accuracy  i  f  a  discontinuity  exists  in  the 
sol ut i on . 

Spectral  — f  i 1  ter i ng ,  solution  smoothing,  and  dissipation 
terms  are  neccessary  to  stabilize  the  pseudo-spectral  scheme 
if  a  discontinuity  is  to  be  "captured". 

Pseudo-spectral  schemes  are  not  suited  -for  the  general  inviscid 
transonic  problem. 

The  proposed  -flow  solver  will  be  uni  versal  to  any  geometry, 
thus  reducing  uncertanties  usually  introduced  by  different 
•flow  salvers  on  di-f-ferent  geometries. 

The  proposed  method  will  also  be  set  up  for  time-accurate 
sol ut ions. 

The  tr i angul ar /tetrahedral  elements  of  the  proposed  method 
will  allow  the  development  of  very  general  contour-plotting, 
streamline  tracing,  and  other  post-processing  graphics  packages. 


REFERENCES 

1.  Lax,  P.D.,  "Hyperbolic  Systems  o-f  Conservation  Laws  and  the 
Mathematical  Theory  o-f  Shock  Waves",  Vol.  11,  Regional  Conference 
Series  in  Applied  Mathematics,  1973,  SIAM,  Phi  1 adel phi  a. 

2.  MacCormack,  R.W. ,  and  Paullay,  A.J.,  "Computational  Efficiency 
Achieved  by  Time  Splitting  of  Finite  Difference  Operators", 

A I A A  Paper  69-354,  Cincinnati,  Ohio,  1969. 

3.  Warming,  R.F.,  and  Beam,  R.M. ,  "Upwind  Second-Order  Difference 
Schemes  and  Applications  in  Aerodynamic  Flows",  AIAA  Journal  , 

Vol.  14. ,  No.  9,  Sep.  1976. 

4.  Cooley,  J.W.,  and  Tukey,  J.W.,  "An  Algorithm  for  the  Machine 
Calculation  of  Complex  Fourier  Series",  Mathematics  of 
Computers,  XIX,  1965,  pp.  45-51. 

5.  Gottlieb,  D. ,  Hussaini,  M. Y. ,  and  Orssag,  S.A. ,  "Theory  and 
Applications  of  Spectral  Methods",  ICASE  Report  83-66,  Dec.  1983. 

6.  Sakel  1  ,  L.  ,  "F'seudospectr al  Solution  of  Inviscid  Flows  with 
Multipule  Discontinuities",  NRL  Memorandum  5147,  Aug.  1983. 

7.  Sakel 1,  L. ,  "Pseudospectr al  Solutions  of  One-  and  Two-Dimensional 
Flows  with  Shock  Waves",  AIAA  Journal ,  Vol.  22,  No.  7,  July  1984. 

3.  Ou,  Y.R.,  Private  Communications,  USC ,  1985-1986. 

9.  Hussaini,  M.Y.,  Kopriva,  D.A.,  Salas,  M.D.,  and  Zang,  T.A., 
"Spectral  Methods  for  the  Euler  Equations:  Part  I  -  Fourier 
Methods  and  Shock-Capturing",  AIAA  Journal,  Vol.  23, 


Jan.  1985,  pp.  64-70. 

10.  Hussaini,  M.Y.,  Kopriva,  D.A.,  Salas,  M.D.,  and  Zang,  T.A., 

"Spectral  Methods  for  the  Euler  Equations:  Part  II  -  Chebyshev 
Methods  and  Shock-Fitting",  AIAA  Journal,  Vol.  23, 

Jan.  1985,  pp.  234-240. 


Tim 

X 

TAtl- 

TAU* 

U* 

SHOCK  SPEED 

•.75999 

9.37599 

9.99999 

•.99999 

9.59999 

9.56999 

I.8HII 

9.39944 

-9. 99956 

9.93671 

9.45349 

9.47679 

•.859ft 

9.42275 

-9.99225 

9.67321 

9.41122 

9.45561 

I.908M 

9.44594 

-9.99496 

9.19991 

9.37237 

9.43618 

9. 95999 

9.46641 

-9.99859 

9.14659 

9.33674 

9.41837 

1.99999 

9.48692 

-9.91398 

9.18278 

9.36414 

9.49297 

1.95999 

9.59666 

-9.91834 

9.21858 

6.27436 

9.38718 

1.19999 

9.52568 

-9.92432 

9.25378 

9.24718 

9.37359 

1.15999 

9.54495 

-9.93995 

9.28828 

9.22249 

9.36129 

1.29999 

9.56183 

-9.93817 

9.32292 

9.19984 

9.34992 

1.25999 

9.57997 

-9.94593 

9.35494 

9.17939 

9.33965 

1.39999 

9.59582 

-9.95418 

9.38699 

9.16963 

9.33932 

1.35999 

9.6121?- 

-9.96288 

9.41817 

9.14367 

9.32183 

1.49999 

9.62892 

-9.97198 

9.44845 

9.12826 

9.31413 

1.45999 

9.64355 

-9.98145 

9.47784 

9.11429 

6.39714 

1.59999 

9.65875 

-9.99125 

9.59633 

9.19161 

9.39981 

1.55999 

9.67365 

-9.19135 

9.53394 

6.99913 

9.29597 

1.69999 

9.68828 

-9.11172 

9.56968 

9.97974 

9.28987 

1.65999 

9.79265 

-9.12235 

9.58658 

9.97935 

9.28517 

1.79999 

9.71689 

-9.13329 

9.61164 

9.96186 

9.28993 

1.75999 

9.73976 

-9.14424 

9.63599 

9.95429 

9.27719 

1.89999 

6.74452 

-9.15548 

9.65937 

9.94731 

9.27365 

1.85999 

6.75813 

-9.16687 

9.68299 

9.94119 

9.27955 

1.99999 

9.77159 

-9.17841 

9.76497 

9.93553 

9.26777 

1.95999 

9.78491 

-9.19999 

9.72534 

9.93955 

9.26527 

2.99999 

9.79812 

-9.29188 

9.74594 

9.92699 

9.26395 

2.95999 

9.81123 

-9.21377 

9.76587 

9.92212 

9.26196 

2.19999 

9.82423 

-9.22577 

9.78518 

9.91869 

6.25939 

2.15999 

9.83716 

-9.23784 

6.89387 

9.91548 

9.25774 

2.29999 

9.85991 

-9.24999 

9.82198 

9.91274 

9.25637 

2.25999 

9.86289 

-9.26229 

9.83953 

9.91934 

9.25517 

2.39999 

9.87553 

-9.27447 

9.85653 

9.99826 

9.25413 

2.35999 

9.88822 

-9.28678 

9.87392 

9.99647 

9.25323 

2.49999 

9.99986 

-9.29914 

9.88991 

9.99494 

9.25247 

2.45999 

9.91347 

-9.31153 

9.99452 

9.99365 

9.25183 

2.59999 

9.92695 

-9.32395 

9.91957 

9.96259 

9.25139 

2.55999 

9.93869 

-9.33649 

9.93418 

9.99173 

9.25987 

2.69999 

9.95113 

-9.34887 

9.94836 

9.99167 

9.25953 

2.65999 

9.96365 

-9.36135 

9.96213 

9.99957 

9.25929 

2.79999 

9.97616 

-9.37384 

9.97552 

9.99924 

9.25912 

2.75999 

9.98867 

-9.38633 

9.98852 

9.99995 

9.2S993 

2.79533 

1.99999 

-9.39766 

1 .••••• 

9.99999 

9.25999 

TABLE  I: 


"Tabulated  Characteristics  o-f  the 
F'ropogating  Shock  during  its 
growth  -from  birth  to  maturity." 
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FIGURE  4:  "Numerical  solution  o-f  Burgers 'equation 

in  Quasi  -Li  near  -form  using 
MacCormack's  method." 
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FIGURE  5: 


'Numerical  solution  o-f  Burgers'  equation 
in  Ful  1  y-Coriservat i ve  -form  using  Beam 
Warming's  implicit  method. 
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Numerical  solution  o*  Burgers  equation 
in  Quasi  -Linear  -form  using^Beam 
Warming's  implicit  method. 
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FIGURE  8:  "Numerical  solution  o-f  Burgers  equation 

in  Quasi -Linear  form  using  Beam-Warming  s 
method  with  dissipation." 
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